Quantitative reliability of Migdal-Eliashberg theory for strong electron-phonon 
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We reassess the validity of Migdal-Eliashberg (ME) theory for coupled electron-phonon systems for 
large couplings A. Although model calculations have found that ME theory breaks down for A ~ 0.5, 
it is routinely applied for A > 1 to strong coupling superconductors. To resolve this discrepancy 
it is important to distinguish between bare parameters, used as input in models, and effective 
parameters, derived from experiments. We show explicitly that ME gives accurate results for the 
critical temperature and the spectral gap for large effective A. This provides quantitative theoretical 
support for the applicability of ME theory to strong coupling conventional superconductors. 

PACS numbers: 74.20.-z,71.10.-w,63.20.Kr 



I. INTRODUCTION 

The theory of conventional superconductivity, where 
pairing is mediated by the couphng of electrons to lattice 
vibrations, is considered as one of the major achievements 
of the twentieth century condensed matter physics. It 
is based on the relatively simple Eliashberg equations^i 
where vertex corrections are neglected. The applica- 
tion of these equations is justified by Migdal's theorem;^ 
which states that vertex corrections are proportional to 
the effective couphng strength A and the ratio of phonon 
cjph and electronic energy scale W, which usually is of 
the order 1/100 and lessi^ii. The properties of the super- 
conducting state in Migdal-Eliashberg (ME) theory are 
then largely determined by A. and a phenomenologi- 
cal parameter for the Coulomb repulsion fi* . ME theory 
is also routinely used as a standard pairing theory for 
other situations where a bosonic pairing mechanism is 
analyzed)^ 

A number of studies mainly based on the Holstein 
model, which go beyond the ME theory, have illustrated 
that vertex corrections cannot be neglected for large cou- 
pling strength even if the ratio ujph/W is very small and 
that the ME theory becomes inaccurate wheni A ex- 
ceeds a certain valuci^"— In the adiabatic limit, Benedetti 
and Zeyher— found a breakdown of Migdal's theorem 
due to the appearance of additional extremal paths in 
the action for A ^ 0.4»i Capone and Ciuchi^i. found 
quantitative deviations of self-consistent ME calculations 
from DMFT already for intermediate coupling strengths 
and qualitatively different behavior for stronger coupling. 
Alexandro\i2, argued that even in the adiabatic limit ME 
theory breaks down due to bipolaron formation and sym- 
metry breaking when A exceeds one. For strong cou- 
pling superconductors values of A of the order 1-3 are 
commonly quotedii^ii The model calculations therefore 
suggest that strong coupling superconductors do not lie 
within the range of applicability of the ME theory. 

The purpose of this paper is to bring the results from 
the model studies in a form that they can be compared 
in a meaningful way to the standard diagrammatic ap- 



proach for superconductivity. This allows us to clarify 
the quantitative reliability of ME theory for relevant val- 
ues of A and ujph/W . If phonon renormalization occurs 
it is necessary to distinguish the bare model parame- 
ters from the effective parameters describing the phonon 
properties ji^"— The latter correspond to the ones derived 
from experiment or density functional calculations. We 
study the Holstein model in the limit where the lattice 
has infinite dimensions. In this limit the dynamical mean 
field theory (DMFT)i^ becomes exact. These DMFT re- 
sults serve as a benchmark for ME calculations. 

We show in qualitative agreement with earlier work 
in the normal phased that self-consistent ME calcula- 
tions become inaccurate already at moderate bare cou- 
pling both for electronic and phonon properties. How- 
ever in contrast with previous interpretations, we show 
that at these bare couplings the effective coupling is very 
large, larger than for strong coupling superconductors. 
For effective couplings relevant for strong coupling su- 
perconductors the ME theory is still accurate. 

In many applications of ME theory the phonons are 
not calculated self-consistently, but taken as an input ei- 
ther from a different calculation or experiment. Then one 
is interested in how accurately electronic properties are 
described by the ME equations for a given phonon spec- 
trum. We can check this explicitely by taking DMFT 
as a benchmark for electronic properties and providing 
the full phonon spectrum as an input for the ME calcu- 
lations (termed ME+ph later). We show that the elec- 
tronic properties are predicted very reliably up to large 
effective coupling strengths within ME+ph calculations, 
i.e. with an accuracy of better than 10%. In this pa- 
per we will neglect the effect of the Coulomb interaction 
usually taken into account via the parameter /i*. 

The paper is structured as follows: In Sec. II, we first 
recall the usual definition of the pairing function a^F{iu) 
and the coupling strength A. Then we give explicit details 
for the DMFT and ME approaches. In Sec. HI, results 
for the comparison of the DMFT and ME calculations 
are shown, followed by the conclusions in Sec. IV. 



2 



II. MODEL AND FORMALISM 

A. Pairing function 

The pairing function a^F{ijj) is an essential ingredi- 
ent for conventional superconductivity. It can be defined 



(1) 



where po is the electronic density of states at the Fermi 
level, gk,k' the electron-phonon coupling matrix element 
and Pq{ijj) the phonon spectral function, related to the 
phonon propagator as 



2uj 



(2) 



These are the dressed phonon quantities of the interact- 
ing system. In conventional theory these are often taken 
from experiment or estimated by a different method, and 
then inserted in the Eliashberg equations to solve for Tc, 
the spectral gap and other properties. As the properties 
of conventional superconductivity are mostly confined to 
a small window around the Fermi energy, often a Fermi 
surface average is used, 

a^F{io) = \J2 <k'F{Lo)5{ek - p)5{£k' - m)- (3) 



Pa 



k,k' 



Then the superconducting state is largely determined 
through the coupling constant^iiiii^ A, 



2 Jduj 





(4) 



B. Holstein model and DMFT approach 

Our quantitative test of the ME theory is based on a 
model, which has been frequently used in the literature, 
the Holstein model, 

H = -^(i»j4,Cj,a+h.c.)+wo^&l&» (5) 

i,j,iT i 



c,- „ creates an electron at lattice site i with spin a, and 



&| a phonon with oscillator frequency luq, ni_u = c\^Ci_cr- 
The electronic density is coupled to an optical phonon 
mode with coupling constant g. We have set the ionic 
mass to 7\/ = 1 in ([5]). The local oscillator displacement 
is related to the bosonic operators by Xi = {hi-^h\)/ \/2cjo, 
where h = 1. 



For the DMFT at T = we solve the effective impurity 
problem with the numerical renormalization groupJ^i^ 
(NRG) adapted to the case with symmetry breaking.— 
For the logarithmic discretization parameter is A = 1.8, 
and we keep about 1000 states at each iteration. The 
initial bosonic Hilbert space is restricted to a maximum 
of 50 states. We use a semi-elliptic density of states 
(DOS) for the electrons po{e) = \fW^^ j (2-k1?) with 
bandwidth W = 4i. In terms of the Hilbert transform 
HT[po](2^) we have for the diagonal Green's function 

Gn(*^„) = AgHT[po](£+) + SGHT[po](e-), (6) 

and for the off-diagonal part 

G2i(iw„) = ApWY\p^\{e+) + Bj.HT[po](£-). (7) 

We have defined Aq = {C2{i^n) + £+{i^n)) / {£+{i^n) ~ 
£_(iw„)), Bg = (C2(«w„) + e_(iiLj„))/(e„(ia;„)-e+(ii:j„)) 
Ap = S2i(«Wn)/(e-i-(«w„) - e_(ia;„)), and Bp = 
S2i(iw„)/(e_(za;„) - e+(iw„)), where 



Ci(iw„) - C2(jw„) 
£± = 7, ± 



(8) 



2'\/(Ci(«'^") + C2(«a;„))2 _ 4E2i(iw„)Si2(iw„), 

with C,i{z) ~ z + p - T,ii{z) and C,2{z) = z - ^ - ^22(2;). 
For the Nambu Green's functions we have Gi2(iw„) = 
G2i(i(jJn) and G22{i^n) = — Gii(— iw„). This implies 
Si2(jw„) = Y.2i{ii^n) and E22(iw„) = -i;ii(-iw„) for 
the self-energies. At half filling Gii(ia;„) and Eii(iw„) 
are imaginary functions, whereas G2i(«w„) and Y,2i{ioJn) 
and D{iujjn) and S]ph(?w,„) are real functions. For the 
semi-elliptic DOS the Hilbert transform is given by 



D 



HT[po](^)- de 



Po(e) 



^(^z-sgn(Im(z))\/l2 



4^2 



(9) 

where the square root of a complex number w is given 
by ^yre^^/^, where ip = [0, 27r), such that the imaginary 
part of y/w is positive. 

At finite temperature, we use the continuous-time 
quantum Monte Carlo (QMC) method developed for 
electron-phonon systems<^ To calculate we study the 
susceptibility in the pairing channel x(9i *w„) in the limit 
qr — >■ and iw„ — >■ 0. It can be expressed in terms of the 
irreducible vertex in the particle-particle channel F'pp^ 
which is calculated in the QMC procedure.— It is then 
sufficient to analyze when the largest eigenvalue of the 
symmetric matrix, 

Mn,,n, = i ) [r^PP' {iuJm , iuJn, ; 0)] («w„, ) , 

(10) 

exceeds one. We have defined the pair propagator 
G(ia;„J - G(-iw„J 



X («w„J 



C(-iw„J - C(«Wt,i) 



(11) 



(^{iuJn) = iuJn + P — S(ia;„), where G(ia;„^) is the local 
lattice Green's function in the normal state. 
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FIG. 1: Diagrams for diagonal and offdiagonal self-energy in 
ME theory involving the full electronic Green's function G 
and anomalous part F and the full phonon propagator D. 



C. ME approach 

The matrix equation for the electronic self-energies in 
ME theory (see Fig. [T|) reads 

E(ia;„) = ^ 5^'r3G(iw„i + iuJn)T3D{iuJm), (12) 
where 

G^iiLOn)-' ^GliiLu^r' ^U^u„) (13) 

with Gij{iujn) = J2k Gij,kiiuJn) and G%{iujn)^^ = ia;„ - 
T^isk — m)- Notice that we use the form valid for large 
coordination number, such that only local lattice Green's 
functions enter. These are calculated with the same semi- 
elliptic DOS as in Eqs. dH) and ^ for the DMFT calcu- 
lation. Thus by ME theory we mean the diagrammatic 
theory, which neglects all vertex corrections to (fT^ . and 
this is compared to the full DMFT results. No further ap- 
proximation such as assuming a constant density of states 
or large bandwidth are made. The latter are common ap- 
proximations in the literature on ME theoryj^ii whose ac- 
curacy can be analyzed in an expansion in Lo-p^/W . This 
is, however, not the subject of the present paper, where 
we focus entirely on the accuracy of the theory when ver- 
tex corrections are neglected. 

The pairing function reads a'^F{uj) — pog^ (uj) for 
the Holstein model. In the non-interacting limit we have 
Pq (uj) ~ S{uj — ujq) — 5{u! + ujq), which in Eq. Q gives 
Ao = /9o2g^/wo purely in terms of bare parameters. This 
quantity was used in model studies and denoted by Ai^ 
However, A as defined in Eq. ^ is given for the inter- 
acting system;^"— Then the phonons are renormalized 
via the self-energy Sph(ii^m), 

D{iuj,n)-^ = D%iLo,^)-^ - I]ph(iw™), (14) 

where D^{iuj„i) ~ 2ajo/[(ii^m)^ ^ i^o]- The lowest order 
contribution to the phonon self-energy is 

2g2 

Sph(ja;m) = ^ G{iuJn)G{iu!m -f «a;„). (15) 

^ n 

As in Ref. [Ill we will call Eqs. HH) and (HH) self- 
consistent ME approximation. In the superconducting 
state an additional contribution from the off-diagonal 
Green's function could be taken into account, which is 
however small and it will be neglected in the following. 



We define the peak of the interacting phonon spectral 
function p^{oj) as the effective phonon scale Wph = Wq. 
There is then a mapping of the bare dimensionless pa- 
rameters Ao,wo/t to the effective parameters X^Lo^/t. A 
exceeds the bare Ao due to the phonon renormalization, 
cjQ '^O' '^^^ increased lattice fluctuations as 

shown in the identity valid at T = 0, 

oo 

wd= Jduj p^(cj) = 2ujo{x^), (16) 



which is generally larger than one. For a sharply peaked 
phonon spectrum the first moment sum rule, 

oo 

J duj ujpoi^) = 2^0, (17) 

— oo 

implies ~ loq/uq. From an estimate for the phonon 
softening due to the lowest order diagram, wJ/wq = 
^/l — aXo, and Eq. Q one can then obtain the result 
A = Ao/(l — aAo)'^'^'^ In three dimensions a = 2, and 
for a semi-elliptic DOS in the limit of large dimensions 
we have a = 8/3. 

In order to calculate the gap at T = we solve equa- 
tion p2p both by introducing spectral functions and an- 
alytic continuation to the real axis and for comparison 
directly on the imaginary axis. D{iuj„i) can be calcu- 
lated self-consistently via Eq. ([T5|) or taken as an input 
from DMFT calculations. The latter type of calculation 
is termed ME+ph. In the ME theory Tc is calculated by 
first finding the local lattice Green's function G(iaj„) in 
the normal phase using 

E(iw„) = -— ^ G{iLO,n + iujn)D{i^m), (18) 
rn 

and then employing Eq. (|10p with 

^(PP)(^L^„,,^L^„,;0) = -g^D{tuj„, - ioj^,). (19) 

III. RESULTS 

Let us first of all establish how the bare and effective 
quantities are related at T = 0. At half filling for fixed 
ujQ = O.lt, we plot A in Fig. [5] (a) and luq/loq in Fig. [2] 
(b) both as function of Aq. We show the results from self- 
consistent ME theory on the real axis (RA) and on the 
imaginary axis (lA) in comparison with the full DMFT- 
NRG result. 

A increases slowly for Ao < 0.3 up to values around one. 
Then it rises more rapidly close to values of Ao where in 
the normal state a metal to bipolaronic (BP) insulator 
transition had been found at Xq ~ 0.464 (shown as a ver- 
tical line) The behavior is qualitatively similar to the 
analytic estimate above, however, as A§>l/a the latter is 
a substantial overestimate and diverges too quickly. The 
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FIG. 2: (Color online) Comparison of the selfconsistent ME 
and DMFT results for renormalized quantities as a function 
of the bare Ao: (a) The effective coupling A as defined by 
Eq. Q, inset wd as given by Eq. (fT6|) . (b) The ratio of 
renormalized and bare phonon frequency ujq/ujq. 
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FIG. 3: (Color online) Comparison of quantities for the super- 
conducting state as a function of A obtained from DMFT and 
ME calculations with phonon input from DMFT (ME+ph): 
(a) The spectral gap Asp; the inset shows the relative devi- 
ation d. (b) The critical temperature for the onset of super- 
conductivity Tc. 



region of most interest for our purpose is A ~ 1—3, typical 
values for strong coupling superconductors. This corre- 
sponds to Aq ~ 0.3 — 0.37 in terms of bare parameters. 
The values for A obtained in the self-consistent ME the- 
ory compare well to the DMFT results for smaller values 
of Ao < 0.3, and then start to overestimate this quan- 
tity slightly. For values of Ao closer to the BP transition 
self-consistent ME underestimates A. We also compare 
the effective phonon frequency which decreases with Ao 
towards zero when Aq approaches Aq . This quantity com- 
pares well to the DMFT result for a considerable range 
of Ao, but starts to deviate for Ao > 0.38 or A > 4. 

We can also calculate electronic properties like the 
quasiparticle weight z or the offdiagonal self-energy 
I]°^(0) which roughly determines the spectral gap at zero 
temperature, Agp ~ zY,°^{G). Then one finds good agree- 
ment for small coupling and moderate deviations between 
DMFT and self-consistent ME theory in the intermediate 
coupling regime, and close to the bipolaronic transition, 
similar to the results for z and which have been ob- 



tained by Ciuchi and Caponeii in the normal state. 

Our main objective is to test the validity of the ME 
theory at strong coupling. Hence, we compare the results 
for the superconducting properties Agp and Tc obtained 
from the ME+ph calculations with the full DMFT re- 
sults. In Fig. [3] (a) we show Agp as extracted from the 
spectral function computed from ME+ph calculations on 
the real axis and the corresponding DMFT result. Notice 
that the results are plotted as a function of A now. 

We find very good agreement for small values of A < 1, 
then a regime where ME+ph slightly underestimates the 
value for the gap, before it exceeds the DMFT result for 
larger values of A. By inspecting the relative deviation 
d = (A,^E _ AD'^''T)/^me plotted as an inset we see 
that there is an agreement in the regime A ~ 1 — 3 better 
than 10%. At very large values of A, Agp from ME+ph 
increases stronger than the DMFT result. For similar 
parameters we have also calculated the critical tempera- 
ture Tc as deduced from the Bethe-Salpeter equation of 
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FIG. 4: (Color online) Guideline for the quantitative reliabil- 
ity of ME theory. The diagram shows points (diamonds) in 
the tJo/W'-A-plane where the deviation d between the DMFT 
and ME+ph theory is d ~ 0.1. The dashed line is given by the 
functional form ujq/W — c/A, which follows from an analyt- 
ical estimate where vertex correction exceed a certain value, 
c = 0.05 was used. 

the uniform pair susceptibility)^ The comparison of ME 
theory and DMFT-QMC result is shown in Fig. [3] (b). 
Good agreement is found in the relevant range for A. 
DMFT-QMC systematically shghtly underestimates the 
phonon renormalization, which accounts partly for the 
discrepancy of Tc from the ME+ph calculations. In a 
related approach Marsiglio found for a 4x4 cluster that 
the self-consistent ME theory agrees well with QMC cal- 
culations for the pairing susceptibility^^ 

By doing similar comparisons for different bare param- 
eters we mapped out for which values of the effective pa- 
rameters A and cjq/W DMFT and ME+ph show good 
agreement, i.e. d < 0.1. The results are shown in Fig. |3] 
and can be well understood in terms of the effective ex- 
pansion parameter of ME theory Xujq/W, which should 



not exceed 0.05 for good accuracy. The results can serve 
as a guideline for the application of ME theory with re- 
liable phonon input. 



IV. CONCLUSIONS 

We have assessed the validity of the ME theory. We 
calculated accurately the effective coupling strength A in 
terms of the bare coupling strength Aq . For intermediate 
Ao the system is close to a bipolaronic metal-insulator 
transition and A is very strongly enhanced. Close to this 
point ME theory breaks down. However, for typical val- 
ues for strong coupling superconductors, A '--^ 1 — 3, the 
ME theory is very accurate for small values of u!p\i/W. 
This result is demonstrated explicitly for the Holstein 
model in the limit of large dimensions, where most of the 
spectral weight of the pairing function is located at Wph. 
We expect that this result is also applicable for more 
general forms of pairing functions a'^F{uj) in three di- 
mensions with an appropriate cut-off scale Wph. In many 
applications of ME theory a momentum average over the 
Fermi surface is taken, such that the situation is similar 
to the one studied here. However, the momentum de- 
pendence can be important in certain cases especially for 
lower dimensional materials. For instance in Ref. |25| . the 
momentum dependence of vertex corrections and their 
effect on Tc was analyzed. 
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